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We analyze the disorder limited motion of quantum vortices in a two-dimensional bosonic super¬ 
fluid with a large healing length. It is shown that the excitations of low-energy degrees of freedom 
associated with the non-analytic reconstruction of the vortex core [Ann. Phys. 346, 195 (2014)] 
determine strong non-linear effects in the vortex transport at velocities much smaller than Landau’s 
critical velocity. Experiments are suggested to verify our predictions. 
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Introduction- Entropy production and dissipation in 
superfluids and superconductors are associated with the 
dynamics of vortices. The cores of these vortices play a 
central role in these processes. For instance, the Bardeen- 
Stephen friction force acting on a vortex in type II su¬ 
perconductors comes essentially from the current flowing 
through its normal core Beyond the linear response 
regime, Larkin and Ovchinnikov showed that the quasi¬ 
particle distribution within the vortex core deviates from 
the equilibrium distribution due to the long inelastic re¬ 
laxation time. The heating in the core of the vortex leads 
to a nonlinear behavior of the friction force acting on the 
vortex even when the value of the applied current is much 
smaller than the critical one mg. 

There is no analogous scenario for vortex flow in 
bosonic superfluids. This is because, unlike superconduc¬ 
tors (and ^He fermionic superfluid), quantum vortices in 
such a superfluid were thought to posses an essentially 
featureless cores, see e.g. Ref. [5]. 

However, we showed recently [g that vortices in two- 
dimensional bosonic superfluids experience non-analytic 
reconstruction of their cores when moving with respect to 
the flow. The theory [g developed for the limit, 1, 

where n is the bosonic density and ^ is the healing length 
(i.e. the size of the vortex core), predicts the following: 
i) The low energy degrees of freedom are not exhausted 
by the position of the vortex itself but must include the 
precession of the vortex around its guiding center; ii) This 
precession can be characterized by the kinetic momenta 
p= {px,Py) such that 


classical quantization Hk(jP) gives discreet energy levels 
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where Z > 0 is an integer; v) Excited states {I > 1) decay 
due to the phonon emission but the relaxation time. Tin, 
is large and discrete levels are distinguishable 
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The spectrum (lb) with vortex mass ~ was 


discussed extensively in the literature concluding that the 
excited states of the vortex are not relevant for the low- 
energy dynamics. The large logarithmic factor in the 
mass, see Eq. makes the dynamics much slower than 
it was previously thought [g. 



[Px,Py] = i2.TT(Tnt? , (la) 

where n is the bosonic density far from the vortex core. 


and cr = ±1 is the vorticity, Eq. (la) corresponds to 


the Lorentz (Magnus) force acting on a moving quantum 
vortex with respect to the superfluid; iii) The momentum 
dependence of the kinetic energy is non-analytic 
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where m is the boson mass, and a = 0.802 ...; iv) Semi- 


FIG. 1. Setups for a) cold atom (BEC) based and b) helium 
films superfluids. In BECs, vortex (black dot) can be intro¬ 
duced by phase imprinting [g, while motion is induced by 
sweeping (with velocity v) a disordered potential generated 
by a speckle pattern (brown impurities) [g. For helium films 
vortices can be introduced by rotating the sample [g (so that 
the vortex lattice is not yet formed), and the superfluid mo¬ 
tion is obtained by evaporative heater located at one end of 
the system and a reservoir of superfluid helium at the other 
end . 

Whereas a direct observation of the core dynamics is 
difficult, its manifestation via the dissipative motion of 
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vortices, like those in superconductors mm , is experi¬ 
mentally accessible, as we demonstrate in this Letter. 

In a stationary clean superfluid the vortices move with 
the flow and the internal degrees of freedom are not ex¬ 
cited. The experiments to pinpoint the excitations in 
disordered systems are sketched in Fig. For a BEC 
[Fig#)] the motion of the vortices can be observed in 
situ whereas for the helium film realization [Fig. [^)] it 
can be deduced from measurements of the chemical po¬ 
tential gradient along the superfluid flow by means of 
differential pressure transducer [^. The superfluid flow 
is obtained by evaporative heater located at one end of 
the system and a reservoir of superfluid helium at the 
other end [lOj . Vortices with the same vorticity are in¬ 
duced by mounting the sample in a rotating cryostat, see 
e.g. Ref. [TT] . 

The main quantities are the non-linear susceptabilities: 


vf = Xf{v)v] ve = cfxe{v)v, ( 2 ) 

where the vortex, ve,f, and superfluid velocities, v, are 
introduced in Fig.[^('(;’s are to be understood as averages 
over the disorder realizations). Our predictions for the 
susceptabilities are summarized in Fig. and the analytic 
expressions are given at the end of this paper. 

In BEC systems, the motion of the vortex can be im¬ 
aged and xf,e are directly measurable. For helium hlms 
the chemical potential gradient, V/r, satishes the relation 

= -2TTnNvXE(v)v, (3) 


where Ny is the density of the vortices per unit area. 
Thus xe{v) can be extracted from the ratio of the chem¬ 
ical potential gradient to the superfluid velocity. More¬ 
over, the N-shape of xe{v)v leads to the situation where 
the same vortex current can be realized for three possible 
value of the supercurrent v, thus leading to the filament 
instability of the supercurrent itself in large samples. 

Qualitative discussion - To discuss the motion of the 
vortex it is necessary to supplement the kinetic energy 


(lb) with the fields coming from the motion of the su¬ 


perfluid surrounding the vortex. The corresponding ef¬ 
fective theory [5] is conveniently written within Popov’s 
formalism [12] mapping the problem of two-dimensional 
superfluid to two-dimensional nonlinear electrodynamics. 
In this mapping, vortices become charged particles with 
charge cr27rfi, the electric field E is related to the super¬ 
fluid current j as E = —ej (where e is the antisymmetric 
tensor of the second rank acting on the spatial coordi¬ 
nates), and the magnetic field, B, is the boson density, 
n. In Popov’s variables, the effective Hamiltonian of a 
vortex (which includes the reconstruction of the core) is 

' ri- peE nh^Btr) 

Ti = +2'Kahip + -1-, (4) 

B m 

where p = V ahA{r,t) is the kinetic momentum 
of the vortex (la), while V is the canonical one, (p{f,t) 



X 


FIG. 2. Non-linear susceptabilities as the function of the 
velocity v expressed in terms of the “threshold” velocity 
V* determined by the disorder and the phonon temperature 
X = v/v*. Different curves correspond to different phonon 
temperature T expressed in units of the disorder energy e^, 
see Eq. (7aI, 2 = Tjtd- Naturally, the threshold behavior 
becomes more pronounced with the lowering of the tempera¬ 
ture. 


and A(f, t) are the scalar and the vector potentials, re¬ 
spectively. The physical fields are E = —Vw — dtA and 
B — V X A. The gauge invariance of Eq. (M is nothing 
but vorticity conservation. 

The third and the fourth terms in Eq. (|^ are not 
present in usual electrodynamics. The third term ex¬ 
presses the fact the vortex executes its motion around 
the guiding center moving together with the superfluid 
and manifests Galilean invariance. The last term is the 
energy of the core depending on the superfluid density 
outside the core, B(r). We will see shortly that this term 
is important for the scattering of the vortex by disorder. 

It is instructive to write the equation of motion from 
the Hamiltonian (|^. Suppressing the non-important ef¬ 
fects of the spatial inhomogeneity of E/B, we have 


p = 2'Kah {^E + Bef^ — irh^WB{E)/m 
f= dpHkip) + eE/B. 


and combining Eqs. (5a)-(^5bL we find 


p = Eujc{p^)ep — Fh?VB{r)/m, 


(5a) 

(5b) 


(5c) 


i. e. the electric field cannot excite the external degrees of 
freedom of the vortex. Such excitations can be induced 
only by scattering on inhomogeneities of the bosonic den¬ 
sity that we will describe now. 

The analogy of the problem with the motion of elec¬ 
tron in magnetic field enables us to use the formalism 
developed for the nonlinear magnetotransport in a two 
dimensional electron gas El- Consider the small Gaus¬ 
sian variations in the boson density n = + Sn, 

{Sn{ri)6n(r2))q = 'ynoG{qra), 


( 6 ) 
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where ■ ■)q denotes the disorder averaging and Fourier 
transform over ri — r 2 , the parameter 7 <C 1 describes 
the disorder strength, Tc is the correlation radius [13] , and 
G{x) is the dimensionless function whose precise form is 
not important for us provided that it drops fast enough at 
X ^ 1. From Eq. (laI it follows that > h?no ^ 
which means that the vortex experiences small angle scat¬ 
tering by the disorder. The relaxation times can be es¬ 
timated by treating the last term in Eq. within the 
Fermi golden rule and neglecting the curving of the tra¬ 
jectories between scattering events. It gives the following 
relaxation time for the momentum direction: 


1/Ttr(e) = UJc {cdlef^'^ , 


(7a) 


where is the characteristic energy scale [TS] associated 
with the disorder. If the kinetic energy of the vortex e 
is larger than the vortex precesses many times before 
changing its position, otherwise the vortex scatters into 
the new position before it manages to complete the circle. 

There is another time scale describing the scattering at 
all angles (and not only those which change the direction 
of the momentum significantly) [15j : 

1/T,(e) = Wc (cq/e)^^^ , Eg » Cd. (7b) 


If uJcTq > 1, one can neglect the interference associ¬ 
ated with the coming back to the same scattering cen¬ 
ter (Shubnikov-de-Haas effect). We will assume e < €q, 
whereas the relation between and e may be arbitrary. 

Consider now a vortex in the coordinate frame moving 
with velocity v = eE/B [moving disorder in this frame 
does not change Hk because of the last term in Eq. @]- 
To be at rest in the laboratory frame (which would be 
consistent with disorder pinning), a vortex should have 
the directed velocity —v. If there were no disorder, Ttr —>■ 
00 , this directed velocity would process and average to 
zero. The presence of disorder allows for rotation by 
small angle uJcTtr ^1. As a result the velocity acquires 
a component along E, ve — cruJcTtrE/B. The motion 
along the electric field leads to Joule heating, and the 
power produced by vortices with typical energy e < Cd 


[see Eq. (7a)] is 


P(e)/(27rS) ~ avEE ~ UcTtrE^lB ~ E^/B (e/cd) ^ 


3/2 


In the opposite limit uJcTtr ^ 1, the circular motion aver¬ 
ages out the dissipative current, and only rare scattering 
events contribute to the dissipation power. Therefore 
the dissipative current should be proportional to 1/rtr, 
and on dimensionality grounds it leads the replacement, 
UJcTtr l/{uJcTtr), 1.6. for e > Crf, 


particular form. If there were no inelastic processes the 
distribution function of the vortices in the energy space 


/(e) would never be stationary. Phonon emissions (Ic) 


remove the energy from the vortex core, and the extra 
energy accumulated by the vortex with reference to the 
starting energy e can be estimated as 


A(e) ~ {2-KhTinE'^/B) (cd /(e^ + e^) 


( 8 ) 


Therefore, for large enough E, there exists a region where 
A(e) > e. The distribution function in such a region is 
almost constant, see Eig. Non-equilibrium currents, 
however, are determined by the energy derivative of the 
distribution function {—df/de), shown in inset of Eig. 
Thus, the currents are not determined by the whole dis¬ 
tribution function but only by regions at small and large 
energies where the dissipative currents are suppressed. 
This explains the non-linearity of dissipation as function 
of E/B, and the drop in the dissipative current at large 
E/B. The quantitative qualitative requires full kinetic 
description of the problem outlined below. 



FIG. 3. The shape of the non-equilibrium distribution func¬ 
tion /(e) calculated from Eq. ( |15[ ). The steep drop at low en¬ 
ergies is the thermal distribution and the high energy tail orig¬ 
inates for the energy resolved “heating” Q. For ei < e < e 2 , 
the relation A(e) > e holds. 

The kinetic equation has the standard form in the 
energy-angle variables m- Suppressing the spatial de¬ 
pendence of the distribution function, /(e, (/>), we obtain 


dtf + aujc{e)d^f = Stci[f] + Sti„[/]. 


(9a) 


The collision integrals in the right-hand-side of Eq. (9a) 
describe probabilistic processes. The disorder generates 
small angle scattering (angular diffusion): 


Stelif] 


’ d 

d E ■ p 

1 

' d 

E ■ p d 

d(j) 

de B 


d(j) 

B de 


f, (9b) 


P{e)/{2TrH) ~ uJcTtrE^B ~ E^B {cd/ef^ , 

Thus the generation of energy is a peaked function 
of e, and non-equilibrium effects are associated with its 


wher e t he elastic relaxation time Ttr{c) is given by 
Eq. (7a). We defined the angle (j) so that E ■ p = 
Ep{e) cos (j), and p(e) = ^2eM„(e), see Eq. (lb). The 
extra term in addition to the angular derivative is the 
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Galilean correction to the vortex energy in the moving 
superfluid [third term in Eq. @]- The energy transfer 
in the phonon emission is small and can be described by 
Focker-Planck terms. Neglecting the effects of the field 
E and the disorder on the inelastic collision, we obtain 


Stinlf] — 


.(e) 


1 +^ 1 - 

oe 


f, 


(9c) 


where T is the phonon temperature, and the inelastic 
rate is given by Eq. @. 

The vortex velocities (in the convention of Fig. are 


.de 


E 


.de 


ve= f ^ cos 4> d4>de; vp= —+ f—sin 4>d(l)de, (10) 


dp 


B 


dp 


with the normalization J d(j)def = 1. 

Solution of the kinetic equation proceeds in a standard 
way e.g. considering the heating effects in metals. Let 
= /o(e) + /i(e;<('), and / dcffiie; (f) = 0. The an¬ 
gular dependent part of the distribution function /i is 
massive and can be found to first order in the perturba¬ 
tion in E. For the same reason, inelastic collision effects 
on fi can be also neglected, and we find 


fi = 


Ep{e) auJcie)Ttr{e) cos (j) — sin (p 


B 


1 + Ulcie^Ttricy 


dfo 
' de 


( 11 ) 


Substituting /(e; p) = /o(e) + /i(e; (p) back into Eq. (9a), 


using Eq. ( |11[ ) and integrating the result over the angle p 
we obtain the spectral diffusion equation dtfo + dj^ = 0, 
where the spectral flow current is given by 


je = -e/o(e)Mn(e) - D,dJo{e). 


(12a) 


The spectral diffusion is caused both by the inelastic pro¬ 
cesses and by Joule heating due to the electric field: 


D, = 


eT 

m(e) 


Epje) 

B 


UJcic)‘^Ttr{e) 

1 + uJcie)^Ttr(e)^' 


(12b) 


In the stationary state the spectral flow is absent = 0 
and we obtain from Eqs. (12) 


dfo fo 
de T^ffie) 


— 0; Teff — T+ 


2TrhTinE'^ {ed ef^^ 


B 


, ( 13 ) 


where we used p^ = 2M^e, Wc = 2TThB/My, and the 
explicit energy depedence of the transport relaxation rate 
(7a). The meaning of the last term in the expression for 
the effective temperature Te// has been already discussed 
in derivation of Eq. ^ . 

Equations (10) and (|Il|) give 


ve 

Vp 


E 

B 


± {ed e) 


3 / 2 ' 


dfo 
' de 


ede 


The normalized solution of Eq. (13) is 
F{e) 


fo{e) = 


Cde.Eie,) 


F = exp 


r dei 

lo di'effiei) 


(14) 


(15) 


Substituting Eq. (15) into Eq. (14), restoring v = E/B, 


and matching overlapping asymptotes for the integrals 
we obtain for the susceptabilites Xv=e,f of Eq. ([^: 


Xv = 


•itIF" di 


(16a) 


where the exponents are <5^ = —2/5; Sf. = 4/5; 5'^ = 
5/2; = 4, and the numerical prefactors are all of 

the order of unity: = 0.58..., dfp = 3/2, dp = 

0.583 ..., dp = 3.53_ The nonlinearity occurs at 

“threshold” velocity 

1/2 

, (16b) 

where fd = 327r^/(81-\/3) = 2.25..., and c = h/{mff) 
is the speed of sound. With the logartithmic accuracy, 
one uses ei ^ T > hwc in expressions ([I^-([Td|). At 
smaller temperatures, one should replace T —)■ huic{ei = 


(3 m 

hujc 

1/4 

'mi'^ed 

3/4 

1 


T 






hiOc)- Each fraction in Eq. (16b) is small so that the 
non-linearity occurs at a superfluid velocity much smaller 
than Landau’s critical value. This value of can be 
understood from the condition T^f f{e) = e = 2T, whose 
meaning is obvious from the qualitative discussion and 

Fig.i 

In conclusion, we constructed the theory for the mo¬ 
tion of quantum vortices in disordered two-dimensional 
bosonic superfluids. The excitations of low energy de¬ 
grees of freedom, associated with core reconstruction [6], 
lead to non-linear transport phenomena, see Eqs. ( |16[ ) 
and Fig. resembling those in superconductors [HE]. 
The confirmation of the peak effect in the dissipation 
and the threshold behaviour in the drift provides evi¬ 
dence for the existence of the vortex core reconstruction 
and further our understanding of the dissipative vortex 
transport. 
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